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ABSTRACT 

Over the last decade, substantial efforts have been devoted to understanding the stability properties, transport 
phenomena, and long-term evolution of weakly collisional, magnetized plasmas which are stratified in temper- 
ature. The insights gained via these studies have led to a significant improvement of our understanding of the 
processes that determine the physical evolution and observational properties of the intracluster medium (ICM) 
permeating galaxy clusters. These studies have been carried out under the assumption that the ICM is a homo- 
geneous medium. This, however, might not be a good approximation if heavy elements are able to sediment in 
the inner region of the galaxy cluster. Motivated by the need to obtain a more complete picture of the dynamical 
properties of the ICM, we analyze the stability of a weakly collisional, magnetized plane-parallel atmosphere 
which is stratified in both temperature and composition. This allows us to discuss for the first time the dy- 
namics of weakly collisional environments where heat conduction, momentum transport, and ion-diffusion are 
anisotropic with respect to the direction of the magnetic field. We show that, depending on the relative signs 
and magnitudes of the gradients in the temperature and the mean molecular weight, the plasma can be subject 
to a wide variety of unstable modes which include modifications to the magnetothermal instability (MTI), the 
heat-flux-driven buoyancy instability (HBI), and overstable gravity modes previously studied in homogeneous 
media. We also find that there are new modes which are driven by heat conduction and particle diffusion. 
We discuss the astrophysical implications of our findings for a representative galaxy cluster where helium has 
sedimented. Our findings suggest that the core insulation that results from the magnetic field configurations 
that arise as a natural consequence of the HBI, which would be MTI stable in a homogeneous medium, could 
be alleviated if the mean molecular weight gradient is steep enough, i.e., (V/i)//i > (VT)/T. This study 
constitutes a first step toward understanding the interaction between magnetic turbulence and the diffusion of 
heavy elements, and its consequences for the long-term evolution and observational signatures of the ICM in 
galaxy clusters. 

Subject headings: galaxies: clusters: intracluster medium — instabilities — magnetohydrodynamics 
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1. INTRODUCTION 

Despite the fact that magnetic fields in galaxy clusters are 
too weak to be mechanically important, they can play a fun- 
damental role in the dynamical stability of the dilute gas 
by channeling the transport of heat, momentum, and parti- 
cles. The weakly collisional character of the hot intraclus- 
ter medium (ICM), which is genetically characterized by sta- 
ble entropy gradients according to Schwarz schild's criterion 
(IPiffaretti et al.H2005b ICavagnolo et al.1120091) . enables the ac- 
tion of ma gnetic instabilities th at are sensitive to temperature 
gradients (Balbus 2000, 2004). In particular, the magneto- 
thermal instability (MTI) exhibits the fastest growing modes 
when magnetic field lines are orthogonal to a temperatur e 
gradient parallel to the gravitational field (IBalbusI 120011) . 
whereas the heat-flux-driven buoyancy instability (HBI) does 
so when magnetic field lines are parallel to a temper ature gra- 
dient which is anti-parallel to the gravitational field (IQuataertl 
12003) . 

While the landscape of thermal instabilities that render 
homog eneous, dilu te plasmas unstable has been well ex- 
plored (Kunz 2 01 II) . and even extended to account for the ef- 
fects of cosmic rays (Chandran & Dennis 2006; Sharma et al. 
120101) . very little is known about the effects that composi- 
tion gradients can have on the stability of the dilute ICM. 
If m agnetic fields do not preven t the efficient diffusion of 
ions (Naravan & Medvedevll200Tt IChuzhov & Nusserl 120031: 
IChuzhoy & Loeb 2004]) then the gradients in mean molecu- 
lar weight can be as important as the gradients in tempera- 



ture (see SectionlTlandlOin & Wu 2000: IPeng & Nagaill2009t 
IShtvkovskiv & Gilfanovll2010t IBulbul et alj 120111) and pro- 
vide another source of free energy to feed instabilities. In 
order to obtain a more complete picture of the stability prop- 
erties of the ICM, it is thus important to relax the assumption 
of a homogeneous medium. 

As a first step toward understanding the role of composition 
gradients in the stability of dilute plasmas, such as the ICM, 
we analyze the stability of a weakly magnetized plane-parallel 
atmosphere where magnetic fields play a key role by channel- 
ing the conduction of heat, transport of momentum, and the 
diffusion of ions. Our a nalysis ge neralizes previo us studies 
on the MTI (IBalbusl2001l). the HBI JOuataertl2008l) . and over- 
stable gravity modes fealbus & Revnoldsll2010l) . and reveals 
the subtle roles played by the temperature and the composi- 
tion gradients in determining the stability of the plasma. 

The outline of the paper is as follows. In Section [2] we de- 
scribe the plasma model for a dilute binary mixture of ions. In 
Section[3]we perform the linear mode analysis and we obtain 
the general dispersion relation that governs the linear dynam- 
ics of a weakly magnetized medium which is stratified in tem- 
perature and composition. We analyze in detail the stability 
of the plasma in the regimes where conduction across a given 
scale is, respectively, fast and slow compared to the dynam- 
ical timescale in Sections |4] and [5] In Section [6] we describe 
the physics driving the most relevant instabilities. We discuss 
the astrophysical implications of this study in Section|7] 
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Figure 1. Schematic representation of the geometry involved in the stabil- 
ity analysis of the dilute, magnetized plane-parallel atmosphere with density, 
temperature, and composition gradients. The symbols || and _L label the di- 
rections parallel and perpendicular to the magnetic field, which is assumed to 
lie on the x—z plane, without loss of generality. 



2. MODEL FOR THE MULTI-COMPONENT, DILUTE 
ATMOSPHERE 

2.1. General Considerations for the Plasma Model 

In order to highlight the physical phenomena that emerge 
when composition gradients are accounted for, we focus our 
attention on a dilute binary mixture (e.g., hydrogen and he- 
lium£] in a fixed gravitational field described by3 



| + v.(H = o, 

— = Vx(vxB), 



pa, 



(i) 

(2) 
(3) 



^ I i ( mpp-) = ( P x-p u )iin^ r -v.g 1 , (4) 



dc 
~dt 



V-Q c 



(5) 



Here, the Lagrangian and Eulerian derivatives are related via 
d/dt = d/dt+v-V, p is the mass density, v is the fluid veloc- 
ity, g is the gravitational acceleration, 7 is the adiabatic index, 
and I stands for the 3x3 identity matrix. The symbols _L and 
|| refer respectively to the directions perpendicular and paral- 
lel to the magnetic field B, see Figure \T\ whose direction is 

given by the versor b = B/B = (b x , 0,b z ). The first term on 
the right-hand side of Equation © accounts for entropy pro- 
duction due to viscous heating in a we akly collisional magne- 
tized plasma (see, e.g.,|Hpllweg 1983). These equations have 
been considered in previous works investigating the dynamics 
of the weakly collisional ICM with a single ion species, i.e., 
in the case whe re the concentration of the other i o n species is 
c = 0, see, e.g.. lKunz! i l2inih : IParTish^ et alj 

(120121) . and references therein. 

Equations (d}-© describe the dynamics of a dilute binary 

1 This analysis can be generalized to consider N species. 

2 For the sake of simplicity, we do not consider here the effects of thermo- 
diffusion, baro-diffusion, etc. ((Landau & Lifshitz 1959ft . 



mixture in the low-collisionality regime and they differ from 
standard MHD in three important respects. 

(i) In a weakly collisional magnetized plasma the pressure 
tensor P = p±\ + (pn — p±_)bb is anisotropic. If the fre- 
quency of ion collisions v„ in single ion species magneto- 
fluid is large compared to the rate of change d/dt of all the 
fields involved, then the anisotropic part of the pressure ten- 
sor is small compared to its isotropic part P = 2p± /3 + p\\ /3 
and (see, e.g.. IHollwedl 19851: ISchekochihin et~aTll2005[) 



jgjj ~P±\ 
P 



< 1 



(6) 



The anisotropic component of the pressure tensor in the mo- 
mentum equation gives rise to the phenomenon known as Bra- 
ginskii viscosity. For small pressure anisotropyfl this contri- 
bution is usually written as 



P\\ ~P-L = 3? ?o ( bb - i| 



(7) 



where ijo is the la rgest of the coeffic ients in the viscous stress 
tensor derived by Braginskii (1965). In order to account for 
the effects of collisions between ions of different species in 
the binary mixture, we replace the vu by an effective ion-ion 
collision frequency v°f , which we define in AppendixlAl 

(h) Heat flows mainly along magnetic field lines, because 
the electron mean free path is large compared to its Larmor 
radius. This process is modeled by the second term on the 
right-hand side of Equation © via 



-XtV)T, 



(8) 



where T is the plasma temperature, assumed to be the same 
for ions and electrons, and y is the thermal conductivity pre - 
dominately due to electrons (Spitzer 1962; Braginskii 19651) . 



X 



6 x 1CT 7 T 5 / 2 erg cm" 1 s" 1 K" 



(9) 



(Hi) The composition of fluid elements can change due to 
particle fluxes. Considering the flux of particles 



-£>GV) 



(10) 



on the right-hand side of Equation (O ensures that the diffu- 
sion of ions is mainly along magnetic field lines. This is a 
good approximation when the plasma is dilute enough for the 
ion mean free path to be large compared to the ion Larmor 
radius. The concentration c is related to the mean molecular 
weight [i via 



l s(1 _ c) (l±^) +c (l±^) 



(11) 



where /i; and Z%, with i = 1,2, are the molecular weights and 
the atomic numbers for the two ion species. The isotropic part 
of the pressure tensor is thus 



P 



(12) 



where fee is the Boltzmann constant and itih is the atomic 
mass unit. 



3 In principle, appropriate closure approximations, see, e.g. JSnvder et al.l 
1 1997) may be adopted to address the high- collisionality regim e starting from 
the Chew-Goldberger-Low, or CGL, limit (Chew et al. 1956). A detailed ap- 
plication of this app roach can be found in the context of accretion disks in 
Sharmaet al. (20®. 
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2.2. Initial Background State 

We consider a weakly magnetized, plane-parallel atmo- 
sphere in a constant gravitational field g = —gz. The back- 
ground magnetic field is weak enough that the mechanical 
equilibrium of the atmosphere, with scaleheight H, is main- 
tained via dP/dz = —gp. We assume that the medium is 
stratified in density, temperature, and composition along the 
vertical z-direction. In the equilibrium state, all the particles 
in the plasma are assumed to be described by a Maxwellian 
distribution with the same temperature, so that pu = p± ini- 
tially. 

In general, the background heat and particle fluxes do not 

vanish, i.e., fo-VT 7^ and 6-Vc 7^ 0, unless the magnetic 
field and the background gradients are orthogonal. The exis- 
tence of a well-defined steady state, i.e., V-Q s = V-Q c = 0, 
demands that the background fluxes should be linear func- 
tions of the distance along the direction of the magnetic field. 
However, even if this condition is not strictly satisfied, the 
dynamics of the modes that we consider is unlikely to be sig- 
nificantly affected if the local dynamical timescale is short 
compare d to the timesca le in which the entire system evolves 
(see also lOuataertll2008l) . 

2.3. Validity of the Braginskii-MHD Approximation 

If the pressure anisotropy grows beyond \p\\ — p±\/P — 
j8 -1 , where the plasma (3 = v 2 h /v A , v t h = (2-P//o) 1//2 is the 
thermal speed, and va = B/(4irp) 1 / 2 is the Alfven speed, the 
Braginskii-MHD approximation embodied in Equations ([T]i- 
becomes ill-posed. This is because the viscous term in- 
troduced in Equation (0 not only fails to damp all the ki- 
netic energy available at the viscous scale but also triggers 
various fast-growing, micro -scale plasma instabilities , such 
as mirror and firehose (see ISchekochi hin et all 120051 120081 
and references therein). The growth rates of these instabil- 
ities are of the order of 7 ~ fc||U t h|p|| — P±\/P an d thus 
they can dominate the plasma dynamics at very small scales if 
b|| — V±-\/P ^ P^ 1 - This poses a challenge in numerical sim- 
ulations addressing the non-linear dynamics of the Braginskii- 
MHD equations since these instabilities grow formally at the 
grid scale and some pro cedure must be d evised in order to 
capture their effects (see IKunz et aTll2012l for a detailed dis- 
cussion). 

One possibility, which would prevent the micro-instabilities 
from operating at once, is to ignore the effects of pressure 
anisotropics and the associated Braginskii viscosity. This 
was the appro ach followed i n the seminal pa pers on the MTI 
dBalbusllioOil) and the HBI dOuatoertll2008l) . which showed 
that both instabilities grow on the dynamical timescale set by 
w<L n — (H/g) 1 / 2 . However, because the timescales involved 
in viscous processes are only a factor of a few longer than 
the dynamical timescales on which both the HBI and the MTI 
operate, accounting for small pressure anisotropies can affect 
the range of wavenumber over w hich these in stabilities oper- 
ate, as well as their growth rates dKunzH20Tll) . Furthermore, 
the timescales involved in processes related to ion-diffusion 
are only a factor of a few larger than those involved in viscous 
processes (see below). Since our aim is to understand the in- 
terplay of the various processes involved in determining the 
stability of a medium stratified in both temperature and com- 
position, we retain the term accounting for Braginskii viscos- 
ity, Equation (0, in the momentum Equation We argue 
next that Equations <fXJ — provide an adequate framework to 
analyze the dynamics of small amplitude perturbations of the 



stratified atmosphere described. 

The equilibrium background state over which we perform 
the stability analysis is such that p^ = p±, and thus there is 
an initial period of time for which the pressure anisotropy will 
remain sma ll enough th a t these plasma-micro instabilities can 
be ignored. IKunz et all (12012ft estimate that the amplitude to 
which the fluctuations in the magnetic field can grow before 
these instabilities set in, and thus Equations (HJ-dUi remain 
self-consistent, is roughly given by SB^/B ~ H/(pX m { p ), 
where A m f p stands for the mean free path between particle 
collisions. We can estimate this value for the ICM as fol- 
lows. The plasma j3 increases from ~ 10 2 in the inner cluster 
regions to ~ 10 4 in the outer parts, while the ratio H / A m f p 
decreases from 10 3 -10 2 in the cluster core to 10 2 -10 in the 
outer region. Therefore, the ratio H/((3A m [ p ) is larger than 
unity in the central regions of a typical galaxy cluster and 
decreases outward to roughly 10 -2 . We thus conclude that, 
for the sake of performing a linear mode analysis, which is 
only formally valid when the fluctuations of all the physical 
variables are small, e.g., 5B/B <C 1, the Braginskii-MHD 
Equations (Q]i-© describes the problem under consideration 
self-consistently. These micro-instabilities are likely to play 
an important role in the subsequent non-linear dynamics, but 
addressing this regime is beyond the scope of this study. 

3. STABILITY ANALYSIS 

3.1. Linearized Equations 

The modes of interest have associated timescales that are 
long compared to the sound crossing time an d it thus suf- 
fices to work in the Boussinesq approximation (Balbusl l2001t 
IOuataertll2008l) . In this limit, the equations for the linear per- 
turbations S ~ e <?t+ik-x becomeQ 



aSv - 



-9- 



z - ikv th 



S_p 
P 



Sp±_ 
P 



15% 

(3 B 



iknv A ' 



SB 
B~ 



aSB = ik\\B5v , 

Sp N 2 . 

a — = Sv z + 

P 9 



7-1 



7-1 , /dlnT , dlnTSB^ 
-i kk- — - — Sb,+o 



Sp: 
■ — 

P 



7 

din /1 

dz 



dz 



Sv z - Dk\ — 
11 P 



dz B 



(13) 
(14) 



(15) 



, ' d In it „ d In u SB 
+iDk- I -^-5b z +b z 



dz dz B 

Here, we have defined the anisotropic viscosity coefficient 



(16) 



1*4 



2v 



eff 



(17) 



where v^f is an effective collision rate for the binary mixture 
(see AppendixlAl. and the thermal diffusion coefficient, 



x£_ 
p 



(18) 



4 The effects of Braginskii viscosity in the thermal evolution of the plasma, 
which appear in Equation l|4j as proportional to (pa — p±) 2 , are of higher 
order and thus they do not contribute to the linear analysis. 
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We have also introduced the Brunt— Vaisala frequency, N 2 , 
which, in a medium stratified in density, temperature, and 
composition, is given by 



N 2 = l±l nPp -f =g ±\ n 
dz 



7 dz 



P~T 



(19) 



Note that, in agreement with the Boussinesq approximation, 
the velocity perturbations satisfy k-Sv = and the fluctua- 
tions in density, temperature, and mean molecular weight are 
related via 



5p ST dfi 



T 



= 0. 



/' 



(20) 



3.2. Relevant Timescales Across a Mode 



Because of the several physical processes that play a role 
in the stability of the dilute atmosphere, the dispersion rela- 
tion corresponding to Equations (fl3l-([Tol is rather involved. 
It is thus useful to understand the hierarchy of the timescales 
involved in the dynamics of a single Fourier mode in order 
to make sensible approximations. The analysis below ap- 
plies to the range of local modes with wavevectors parallel to 
the magnetic field for which the fluid approach is valid, i.e., 
H- 1 < fci, < AzLor 



< fe|| \J KnfpH < V K n 1 

where we have defined the Knudsen number 

H 



(21) 



(22) 



The inverse timescales characterizing the diffusion of heat, 
momentum, and particles along magnetic field lines are 



i 



(fc-b) 2 3z/|| 
{k-b) 2 D. 



(23) 

(24) 
(25) 

For a given mode, the ratio between these timescales is in- 
dependent of the direction of the wavevector characterizing 
the perturbation and the background magnetic field and is set 
by plasma processes. Because heat conduction is mostly due 
to electrons, while viscous processes are dominated by the 
dynamics of ions, it could be expected that the associated 
timescales would satisfy t" 1 . However, this is not the 

case and a simple estimate leads to t~ ~ 6T" 1 dKunzl201 II) . 
It could be argued that the timescales involved in viscous and 
diffusion processes should be of the same order because it is 
mainly the ion dynamics what determines both of them. A de- 
tailed analysis of the diffusion coefficient for a binary mixture 
of ions (see AppendixlAl shows that t^ 1 ~ 9t~ x for primor- 
dial composition (c ~ 0.25 or fj, ~ 0.6) and decreases toward 
tT 1 ~ St^ 1 for the compositions expected at the inner core 
of galaxy cl usters according to recent models for helium sed- 
imentation (B ulbul et alJI201 ll) . Since we will be mostly con- 
cerned with the two regimes t" 1 S> Wdyn or ^dyn ^ T c 1 ' 
as long as the ratio t^ 1 /t^ 1 is not too small, its particular 
value will not affect our main conclusions, and we will thus 
consider that 1 ~ t~ 1 . 

On the other hand, whether the timescales set by plasma 
processes are fast or slow compared to the dynamical 



timescale = (H/g) 1 ^ 2 depends not only on the wave- 
length of the mode but also on the direction of the wavevector 
characterizing the perturbation with respect t o the backgr ound 
magnetic field. In particular, as shown in (Kunz 20 fj]), the 
timescale characterizing conduction across a mode with par- 
allel wavenumber kn is related to the dynamical timescale via 



lOfcn X m fpH CJdyn , 



(26) 



where we have assumed 7 = 5/3 in Equation (1231 . Thus con- 
duction is faster than the dynamical time, i.e., t~ /k>dyn > 1- 
if fcii (AmfpiT) 1 / 2 > 1/3. If the wavelength of the mode 
is shorter than this by a factor of t~ 1 /t^ 1 ~ 6, e.g., 
k\\ (Amfp-ff) 1 / 2 3> 1, then viscous and diffusive processes 
are also faster than the dynamical timescale. Therefore, as a 
useful approximate criterion, whether fcii (A m f p -ff ) : / 2 is much 
larger or smaller than unity defines whether the timescales 
associated with plasma processes, for that given mode, are 
shorter or longer than the dynamical time. We will thus con- 
sider two different regimes which we refer to as the "fast" and 
"slow" conduction limit, where the timescales associated with 
the modes considered satisfy, respectively, 



r- 1 >r" 1 



r d " 1 » u dyn if fey » (A mfp i7)- 1/2 , (27) 



Wdyn > T c 1 > T v 1 



r d L if fc|, « (X mip Hy^ 2 .(28) 



3.3. The Weak Magnetic Field Limit 

All the timescales related to plasma processes discussed 
above depend only on the direction of a given wavevector with 
respect to the magnetic field. The only time scale that depends 
explicitly on the strength of the field is the one associated with 
the Alfven frequency oja = k-VA- In order to keep the prob- 
lem tractable, and given that we are already dealing with four 
different timescales, we will focus on the case where the mag- 
netic field is so weak that its only physical role is to channel 
the flux of heat and ions. The advantage of this limit is that it 
allows us to address the anisotropic dynamics of the weakly 
collisional magnetized medium without introducing explicitly 
the timescale associated with cja- 

In what follows we focus our attention on modes for 
which magnetic tension is unimportant and thus oja <C 
minjrjT 1 , Wdyn}- This approximation will be valid for 
two different ranges of parallel wavenumbers depending on 
whether fcy (Amfpff) 1 / 2 is much larger or smaller than unity. 
For the modes for which conduction is faster than the dy- 
namical timescale, i.e., fcii (Amfpi/) 1 / 2 3> 1, we must require 
wa <C ^dyn <C t" 1 . Using the definitions Vth = (9H) 1 / 2 , 
fj = v 2 h /v\, and K n = X m f p /H, we obtain that magnetic 
tension is negligible provided that 



and thus 



1 < fc|| \/X mtp H < y/PK n , 



for t c 1 > w dyn if @K n > 1 



(29) 



(30) 



For the modes for which conduction is slow compared to the 
dynamical timescale, i.e., fen (Amfp-ff) 1 ' 2 <C 1, we must re- 



quire u>a <t c < ^dyn- This is satisfied if 
1 1 



lOy/Wr 



(31) 
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and therefore 

lu a ~ for Lu dyn > r" 1 if /3if n > 10~ 2 . (32) 

The plasma (3 ranges from 10 4 in the outskirts of the 
ICM down to 10 2 in the centers of cool core clusters 
dCariHi & Tavlorl 120021) . and the product (3K n ranges from 
10 3 in the outskirts of the ICM decreasing to 10 _1 in the inner 
regions. Thus the effects of magnetic tension can be impor- 
tant in the inner cluster regions. We address the implications 
of neglecting magnetic tension when analyzing the stability 
of the ICM in further detail in Section [7J 

3.4. General Dispersion Relation for the Dilute, Weakly 
Magnetized Medium 

The dispersion relation corres pon ding to the set of equa- 
tions for the linear perturbations (U3t -(ll6t is 



E" 

i=0 



where the coefficients Aj are given by 

A = l, 
A x = t-\ 



A 2 = 



i k l + k D N 2 +T -i T -i 

U2 d c ' 



, [dhiTK dhiuikl + kl) 
A3 = t~ g ' 



dz k 2 dz k 2 
-1 i k l + k y\ N 2 



Aa=t- 1 t- 1 N 2 — 



while the B; read 



B - kl 

k 2 

R — ^-1 -L 



h 2 k 2 

S,EET7V-^+iV 2 ^ 



k 2 

_ -1 -1 
J 3 = T d T c p 



B^^N^+r^N 2 )^ 



k 2 ' 



r — t- 1 ^- 1 at 2 x y 



Here we have defined 



-6 2 fc 2 + fc 2 + 5 2 fc 2 . 



and the two quantities 



7V 2 M =^-ln(I>), 



d . (T 



(33) 



(34) 
(35) 

(36) 



(37) 
(38) 

(39) 
(40) 
(41) 
(42) 
(43) 



K = {l-2b 2 z ){kl + k 2 y ) + 2b x b z k x k Zl (44) 
= b 2 x k 2 -k\ + b 2 x k 2 yl (45) 

(46) 



(47) 
(48) 



which appear naturally when thermal and composition gradi- 
ents are considered. The dispersio n relation (B3l is identical 
to the one derived in (Kunz 2011) in the limit in which 10 a, 
dp/dz, and D vanish. Note that in the limit of a constant 
composition gradient, i.e., dp/dz — > 0, both and N^,^ 

— > gdhiT/dz, which is the logarithmic gradient that plays 
an important role in the stability of a homogeneous, dilute, 
weakly magnetized medium. 

4. THE FAST CONDUCTION LIMIT 

We first consider the stability of the modes for which con- 
duction is faster than the dynamical time, i.e., t" 1 3> LJdyn- 
This is the regime that corresponds to the well-studied HBI 
and MTI. 



4. 1 . Limit of No Ion-Diffusion 

As a first step toward understanding the effects of composi- 
tion gradients in the behavior of the HBI and the MTI, we 
neglect the diffusion of ions along magnetic field lines by 
setting D = 0. Because we are considering the timescales 
for ion-diffusion and viscous processes to be of the same or- 
der, i.e., r" 1 ~ t7 , we also ignore here the effects of vis- 
cosity and set = for consistency. For the modes for 
which t^ 1 3> uidyn, the dispersion relation (l33l yields a (fast) 
decaying solution, a w — r" 1 , together with the two slow 
modes 



-9 



d In T JC dhiji (k 2 



kl) 



dz k 2 



dz 



(49) 



In a homogeneous plasma, these slow modes contain the well- 
known HBI and MTI, depending on the direction of the back- 
ground magnetic field, i.e., 



dlnTk 2 



for b z = 1 , 



dz k 2 



'MTI ■ 



-9- 



d\aTk 2 x + k 2 



dz 



k 2 



(50) 



(51) 



The conditions for the excitation of the HBI and the MTI are 
thus 



dlnT 

dz 
dlnT 

dz 



> HBI-unstable , 
< MTI-unstable . 



(52) 
(53) 



Note that in both cases, the fastest growing modes are those 
with wavevectors perpendicular to the gravitational field. 

4.1.1. Heat- and Particle-Flux Driven Buoyant Instability 
(D = 0) 



Consider a magnetic field parallel to the gravitational field, 
i.e., b z = 1, and thus k 2 + k 2 = k\. The modified version 
of the modes that become HBI-unstable in a homogeneous 
medium is given by Equation d49| ). In the medium stratified 
in composition, these modes become 



■'<-)- 



_dln(T/x) k\ 

dz k 2 



(54) 



Therefore, neglecting viscous and diffusion processes in a 
medium which is stratified in the mean molecular weight 
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Figure 2. Graphic representation of the stability of modes for which conduction is fast compared to the dynamical timescale, i.e., r c S> i^dyn- The various 
panels show the unstable regions (gray) for each of the modes that can be excited when the background magnetic field is parallel (a, b, and c) or perpendicular 
(a!) to the background thermal and composition gradients. The horizontal dotted line in panel (a) represents dhiT / d\n P = (7 — l)/27; and the dashed lines 
correspond dlnT/dlnP = ±d In fi/d In P. Panel (a) shows the region of parameter space which is unstable to the heat- and particle-flux-driven buoyancy 
instability (HPBI), together with the regions that are overstable to gravity modes. If ions can diffuse efficiently along magnetic field lines, i.e., D ^ 0, there are 

unstable modes that can be driven by diffusion, whether Wdyn > t^ 1 — t^ 1 (b) or tT ~ 7V~ > W(jyn (c). Panel (d) shows the region that is unstable to the 
magneto-thermo-compositional instability (MTCI); in this case, the criterion for instability is insensitive to the value of the diffusion coefficient D. 



leads to modes that are unstable to a Heat- and Particle-flux- 
driven Buoyant Instability (HPBI) if 



d In T d\nn 

> -- 



dz 



dz 



HPBI-unstable . 



(55) 



The threshold temperature gradient for instability can be neg- 
ative if the mean molecular weight increases with height. 

The combination of temperature and composition gradients 
that is HPBI-unstable is shown in panel (a) of Figure|2] where, 
for the sake of convenience, we have defined dimensionless 
variables (dln/i/dlnP, d\x\T/d\x\P) in terms of the loga- 
rithmic pressure gradient d In Pj dz = —1/H. 

4. 1 .2. Overstable Modes (D = 0) 



IBalbu s & Reynolds: (120101) showed that temperature gradi- 
ents that are stable to the HBI can nevertheless be subject 
to overstable gravity modes. This result can be extended to 
include non-vanishing composition gradients, i.e., there is a 
range of modes with t~ 1 > Wdyn and a ~ uJdyn that can be- 
come overstable when the heat- and particle-flux-driven buoy- 
ancy instability (HPBI) does not operate. Calculating these 
modes requires retaining higher order terms in the dispersion 
relation, which becomes 



-v 2 



0. 



(56) 



In this regime, we can treat the first and third terms on the 
left-hand side as perturbations and extend the solutions 
to contain corrections of order /r" 1 : 



±i 



h 2 

-N 2 



N 



2rr. 



(57) 



The combination of temperature and composition gradients 
that is subject to overstability is shown in dark gray in panel 
(a) of Figure [2] 

4.1.3. Magneto-Thermo-Compositional Instability (D = 0) 

In order to understand how the MTI is modified in the pres- 
ence of composition gradients we consider a horizontal mag- 
netic field along the x-axis, i.e., b x = 1 and focus on the 
modes for which r^ 1 ^> Wd yn - In the presence of a gradient 
in the mean molecular weight, Equation ( |49] i gives the gener- 
alization of the modes that become MTI-unstable 



-9- 



dln(T/ M ) kl + kl 
dz k 2 



(59) 



Thus a non-vanishing gradient in the mean molecular weight 
sets an upper bound for the temperature gradients that are 
magneto-thermo-compositional instability (MTCI)-unstable 



d In T d In a 

< 



dz 



dz 



MTCI-unstable . 



(60) 



Panel (d) in Figure |2] shows a graphical representation of the 
region of parameter space that is subject to unstable MTCI 
modes in a medium that is stratified in composition and tem- 
perature. 

4.2. Ion-Diffusion Along Magnetic Field Lines 

We now analyze the effects of including ion-diffusion in- 
duced by the background composition gradients. Since, for 
a given mode, viscous and diffusion timescales are of the 



same order, i.e., t 



we also consider the effects of 



anisotropic viscosity for consistency. 



Therefore, modes that are stable according to the HPBI- 
stability criterion, i.e., 7V|. < (Equation l55l ). can become 



overstable if iVy, 



-N 2 . 



< — TV*. In terms of the dimensionless 
variables introduced earlier, these requirements become 



dlnT 
d\nP 



> max 



d In [i 7 — 1 



4.2. 1 . Heat-and Particle-Flux-Driven Buoyant Instability 

Starting from the general dispersion relation in the case 
where b z = 1, it can be seen that, if we consider modes for 
which conduction is faster than any other timescale, there is a 



HPBI-overstable .(58) fast decaying solution a 



which can be used to self- 
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consistently obtain three more modes satisfying 

k 2 \ k 2 

+ £2 Vd T v 

- T d lN T/S = - (61) 

However, the leading order solutions to this dispersion rela- 
tion depend on whether the dynamical timescale is fast or 
slow with respect to the diffusion and viscous timescales 
across the mode. We thus need to consider these two cases 
separately. 

Slow diffusion. The modes for which cJdyn > Ty 1 — t^ 1 
contain the generalization of the HBI for non-vanishing com- 
position gradient^ 



AT 2 ZL 

f 1 fc2 

r2 



(62) 



which grows dynamically if iVf, > 0, or dT/dz > —dfx/dz. 
This is the same condition for the onset of the HPBI in Equa- 
tion d55l l in the absence of diffusion. It is worth mention- 
ing that even if the combination of temperature and compo- 
sition gradients is such that TVf. < (and thus the plasma 
is HPBI-stable), overstable modes can be excited just as in 
the non-diffusive case. In fact, the condition for overstabil- 
ity can be shown to be exactly similar to d58l l. provided that 

Wdyn > max{r c rV" 1 , t~ 1 t~[ }. There also exists a slower 
mode driven by ion-diffusion 



N 2 , 

'd N 2 



(63) 



which grows in the region of parameter space where 
Nt/h/Nth < or, in terms of the background gradients, 

wherever d\xiT/dh\P < |<f ln/x/c£lnP|, as it is shown in 
panel (b) of Figure [2] Note that these modes can grow on 
a diffusion timescale even if the system is HPBI stable, i.e., 
< 0, provided that iV|y M > 0. As we show below, this 
latter requirement becomes the deciding one for those modes 
for which diffusion is not slow compared to the dynamical 
time. 

Fast diffusion. For the modes for which t~ x ~ t^ 1 > 
Wdyn, there are two solutions that decay on the diffusion and 
viscous timescales, i.e., a 
and a third one 



-Tj 1 and a 



- x kl/k 2 , 



N 2 , 



(64) 



which can become unstable if N%, > 0. This region of 

parameter space in temperature and composition gradients is 
unstable to a diffusive version of the Heat- and Particle-flux- 
driven Buoyant Instability (D-HPBI) 



> ——!- D-HPBI-unstable , 



(65) 



dz dz 

and it is shown in panel (c) of Figure [2] 

In summary, we conclude that when finite viscous and dif- 
fusion timescales are considered, there can be unstable modes 

5 When a finite diffusion timescale is considered, i.e., D^O, the term that 
contributes to the generalization o f the HBI is the second term inside brackets 
in the right hand side of Equation 1161 . 



driven by diffusion (with u>dyn > t^ 1 ) even if the temper- 
ature and the composition gradients do not satisfy inequal- 
ity d55l ), i.e., they are stable to the HPBI in the absence of 
diffusion, provided that d In T > d In /x. Note that this is 
the very requirement for the existence of unstable modified 
HBI modes when finite diffusion timescales are relevant (with 



t7 > Wdyn) and differs from the condition (|55l l. We shall 
discuss the physical reason behind this change in condition of 
instability in Section[6] 

4.2.2. Magneto-Thermo-Compositional Instability ffl/OJ 

In order to understand how ion-diffusion driven by a com- 
position gradient affects the MTI, we consider D ^ and 
modes for which r^ 1 S> uJdyn when b x = 1. In this case, 
the dispersion relation factorizes and leads to a fast decaying 
solution a = — r^ 1 together with 



k 2 



For the modes satisfying o;d yn > t v 1 , the three solutions to 
this cubic equation correspond to a decaying viscous mode 



k 2 



1,2 



U2 

s 



(67) 



and a pair of roots that contain the generalization of the MTI 
in the presence of a mean molecular weight gradient and D ^ 




k 2 



(68) 



which grows dynamically if N? 



T/n 



< 0, or dT/dz > d^/dz 



This is the same condition for the onset of the MTCI in 
Equation ( f60b in the absence of diffusion. As we found be- 
fore, a non-vanishing gradient in the mean molecular weight 
sets a upper bound for the temperature gradients that are MTI- 
unstable. This temperature gradient can be positive when the 
composition gradient is negative (i.e., mean molecular weight 
increasing with height). Note that this statement is indepen- 
dent of the value of the diffusion coefficient D, i.e., and thus 
of whether or not ions diffuse effectively along magnetic field 
lines on a dynamical timescale. 

For the modes such that t^ 1 > Wdyn> the three solutions to 
Equation (l66i l correspond to another viscously damped mode 



;.2 

k 2 ' 



(69) 



together with a generalizatio n of the modes identified as 
Alfvenic-MTI in lKunz! ( f20TTh . i.e., 



N. 



T/n 



2T7 1 k± 



±i 



N 2 , -2- 

T /»k ± 



(70) 



These modes grow dynamically if iV|y^ < 0, which corre- 
sponds again to the combination of temperature and composi- 
tion gradients satisfying inequality (f60b - 

5. THE SLOW CONDUCTION LIMIT 
We now analyze the stability of the modes whose associ- 



1 > t- 1 



For 



ated timescales satisfy Wdyn r c 
a homogeneous plasma, these modes encompass the over- 
stable g-modes studied in Ba lbus & Reyn olds (2010). The set 
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Figure 3. Graphic representation of the stability of modes for which conduction is slow compared to the dynamical timescale, i.e., ui^yn ^ T c ■ The 
various panels show the unstable regions (gray) for each of the modes that can be excited when the background magnetic field is parallel (a, b, and c) or 
perpendicular (d and e) to the background thermal and composition gradients. The solid line corresponds to TV 2 = 0; the horizontal dotted line represents 
d InT/d InP = (7 — l)/2j; and the dashed lines correspond to d In T/d In P = ±d In fj,/d In P. For B = B z , gravity modes (a) can become either stable or 
overstable, while modes driven by conduction can become unstable (b). If ions can diffuse efficiently along magnetic field lines, a new type of mode can become 
unstable (c). For B = B x , both gravity modes id) and conduction modes can become unstable (e); while ion-diffusion only leads to decaying modes. 



of Equations (fi3ll-(l20li allows us to address the behavior of 
these, as well as other new modes, in the presence of a non- 
vanishing gradient in the mean molecular weight and account 
self-consistently for the diffusion of ions along magnetic field 
lines. 



5.1. Heat and Ion Diffusion Along Vertical Magnetic Fields 

In the case where the background magnetic field is parallel 
to the gravitational field, b z = 1, the dispersion relation (l33l 
reduces to 

g + aia 3 + a^rr 2 + 0317 + 04 = , (71) 
with coefficients 

P 

1 1 rv I 1 1 

ai=r c +r d +jjT v wr c , 

k 2 f k 2 \ k 2 



k 2 

-1 -L ¥ 2 



k 2 



(73) 
(74) 



The Routh-Hurwitz stability criteria that predict exclu- 
sively negative real parts for the roots of the quartic polyno- 
mial with real coefficients require a\ > 0, a\a,2 — 03 > 0, 
\ > 0, and 04 > 0. The first of these 



aia2d 3 



a 2 a4. 



conditions is trivially satisfied, while the other three imply, 
respectively, 



N 



7V 2 /At <0. 



(75) 
(76) 
(77) 



In the absence of diffusion, only the conditions d75l l and d76*| l 
need to be met in order to ensure stability, while the condition 
d77| i should also be required for finite diffusion timescales. 

To leading order, two of the solutions of Equation (1711 are 



1 /2 

given by a m ±ia 2 + (03 — a\ai) j1a,2, i.e., 



N 2 



(78) 



which correspond to gravity modes. In the absence of a gra- 
dient in the mean molecular weight, these re duce to the g- 
modes discussed in lBalbus & Reyno lds (2010). The third root 
is given by a w —03/02, i.e., 



N 2 

__ 1 Tfj, 

c TV 2 



(79) 



and corresponds to a mode driven by conduction. Assuming 
that A^ 2 > 0, g-modes are overstable if the condition (l76l is 
not satisfied, while conduction modes are unstable if Equa- 
tion d75l l is not fulfilled. The fourth solution consists of a 
mode driven by ion-diffusion 



N 2 , 
N 2 



(80) 



which is unstable if either criterion (fTBl or (l77t is unfulfilled. 

It is useful to understand what types of modes can be ex- 
cited in the different regions of the parameter space spanned 
by the gradients in temperature and composition. The Routh- 
Hurwitz stability criteria take simple forms when expressed 
in terms of the dimensionless gradients defined in terms of 
the pressure. The classical requirement for stability against 
buoyancy, i.e., N 2 > becomes 



d In T d In n 

< 



7' 



1 



dlnP dlnP 7 ' 

while the conditions d75l l and dTol i become, respectively, 

d In T d In u 
>-- 



dlnP 
dlnT 

dlnP 



< 



dlnP 

7-1 
27 



(81) 



(82) 



(83) 



If ions can diffuse along magnetic field lines, in addition to re- 
quiring that the gradients in temperature, pressure, and mean 
molecular weight satisfy the inequalities ( 1821 and (|83l , the 
inequality ( |77| ) must also be satisfied, i.e., 



d In T d In fi 
d\nP > d\nP 



(84) 
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Panel (a) in Figure [3] shows the regions of parameter space 
where g-modes in Equation ( 1781 are stable, overstable, or un- 
stabl e (g ray area). Panel (b) shows that the modes in Equa- 
tion d79l , which are driven by conduction, can be either stable 
or unstable. Note that in the region of parameter space where 
both gravity and conduction modes are overstable/unstable 
they both grow with comparable rates. Panel (c) in Figure [3] 
shows that the modes in Equation d80b , which are driven by 
diffusion, can be either stable or unstable (gray). Their growth 
rates are estimated to be an order of magnitude smaller than 
either g-modes or conduction modes. The importance of these 
diffusion modes resides in that they can become unstable in 
regions of parameter space which are stable against g-modes 
and conduction modes. 

5.2. Heat and Ion Diffusion Along Horizontal Magnetic 

Fields 

If the background magnetic field is perpendicular to the 
thermal and the composition gradients, i.e., b x = 1, the dis- 
persion relation (l33l becomes Q 



cr 4 + &i(7 J 

where the coefficients 

&1~T- 1 + 



&2C + b^a + 64 = . 



k 2 

-1 ft _L 



k y 



k 2 



k 2 

c v k 2 



h 7 



64' 



k 2 

■T- 1 T- 1 M 2 V 



k 2 

k 2 



(85) 

(86) 
(87) 
(88) 
(89) 



are subject to the same considerations employed in deriving 
the approximate expressions for the coefficients a^. 

The Routh-Hurwitz stability criteria require b\ > 0, 61^2 — 
63 > 0, &1&2&3 — — &| >> and 64 > 0. The first condition 
is trivially satisfied, while, in the limit under consideration, 



I.e., W dyn > T c 

respectively, 



> t v 1 , the other three conditions become, 



N 2 

r2 



N* /It >0. 



(90) 
(91) 
(92) 



The inequality d90l l is always satisfied, since it can be written 

as 



7' 



1 



iPp 



dP 

dz 



> 0. 



(93) 



Therefore, the only independent condition required for stabil- 
ity is N T , > 0, or din T/dz > din fi/dz. 

Two of the approximate solutions to the dispersion relation 

1/2 

are given by a ~ ±ib 2 + (03 — &ifr2)/2&2, i-e-, 




N 2 



(94) 



1 In the absence of diffusion, the only result that is modified in this section 



is that the root a 



-T, for D ^ becomes u = for D = 0. 



These correspond to gravity modes, which cannot become 
overstable on account of Equation (|93l l. The other two so- 
lutions correspond to a conduction and a viscous (decaying) 
mode, which are, respectively, 



N 2 
-1 T/n 

N 2 



7.2 



£-2 1 1.2 



(95) 



We conclude that when the magnetic field is perpendicular 
to the temperature and the composition gradients, the stability 
of g-modes requires only that A^ 2 > 0, whereas the stability 
of conduction modes requires also 



d In T dliifi 
d\nP < d\nP 



(96) 



This is illustrated in panel (d) in Figure [3] which shows the 
regions of parameter space where g-modes in Equation d94l 
are stable or unstable (gray). Panel (e) shows that the modes 
in Equation d95l ), which are driven by conduction, can be ei- 
ther stable or unstable. These regions are significantly dif- 
ferent from the corresponding regions in panel (b), for which 
the direction of the background magnetic field is parallel to 
the direction of the temperature and the composition gradi- 
ents. Note that, unlike the case where b z = 1, g-modes and 
conduction modes cannot be simultaneously unstable when 
b x = 1. 

6. THE PHYSICS DRIVING UNSTABLE MODES 

In previous sections we have seen that the criteria for in- 
stability result from an intricate interplay between the ther- 
mal and the composition gradients. The final expressions for 
the inequalities that must be satisfied for the onset of unstable 
modes depend implicitly on whether conduction across the as- 
sociated scale is faster or slower than the dynamical timescale, 
the direction of the background field, and the ability of the 
ions to diffuse along magnetic field lines. We now analyze in 
detail the eigenmodes corresponding to some of the most rel- 
evant instabilities and shed light on the physical phenomena 
that play a role in determining the stability of the plas ma. This 
exerci se is similar in spir it to the ones presented in IQuataertl 
J20081) and lBalbusI J200T1) to highlight the physics of the HBI 
and the MTI but emphasizes the effects of composition gradi- 
ents and ion-diffusion along magnetic field lines. For the sake 
of simplicity we assume that the perturbations under consid- 
eration correspond to modes with k y = 0. In this case, the 
components of the Lagrangian displacement, d£/dt = Sv, 
are related via £ x = -(k z /k x )£ z . 

6.1. Heat- and Particle-Flux Driven Buoyancy Instability 

Let us first consider a background magnetic field with b z = 
1 and focus on perturbations with wavelengths such that the 
associated timescales satisfy r" 1 ^> Wdyn 3> t^ 1 ~ t^ 1 . 
For these modes, ion-diffusion along the magnetic field lines 
is inefficient. In order to understand the effect of a composi- 
tion gradient, we retain the dominant terms in Equations (TT3T > 
and ( |16*| ), which leads to 



dz 



d/i 

dz 



(97) 



Because the fluctuations in density, temperature, and mean 
molecular weight are related via Equation (l20t . this implies 
that the relative change in density of a fluid element which is 
vertically displaced by £ z is given by 



Sp 



d\i\T d In fj. 



dz 



dz 



HPBI. 



(98) 
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Figure 4. Schematic representation of various unstable modes in a weakly magnetized plasma with temperature and composition gradients parallel to the 
gravitational field g = —gz. The three sets of panels (a,b); (c,d); and (e,f) show the modes that are unstable to the heat- and particle-flux-driven buoyancy 
instability (HPBI), the diffusive-HPBI, and the magneto-thermo-compositional instability (MTCI), respectively. The arrows represent the Lagrangian displace- 
ments, assumed to be of the form £ = £ cos(k x x + k z z), with k x = k z . The continuous lines represent the magnetic field lines, which are assumed to be 
parallel (HPBI, D-HPBI) or perpendicular (MTCI) to the gravitational field in the equilibrium state. The gray-scale contours show the temperature and the mean 
molecular weight fluctuations relative to the background gradients, which are shown with arrows indicating their directions in each of the cases considered. 



It is easier to understand the physics behind this equation by 
analyzing first the effects of each of the two terms on the right- 
hand side separately. 

In a homogeneous medium, dfx/dz = 0, the term propor- 
tional to the temperature gradient is responsible for the onset 
of the HBI when dT/dz > 0. This is illustrated in panel (a) 
in Figure|4] which shows that a fluid element that is displaced 
with £ z > (£ z < 0) is effectively heated (cooled) by the in- 
creased (decrease) flux of heat due to the convergence (diver- 
gence) of field lines that results from the displacement of the 
fluid element. This causes the fluid element to expand (con- 
tract), and thus attain a density which is lower (higher) than 
the surrounding medium. This leads to the runaway process 
known as the HBI. 

An isothermal environment, dT/dz = 0, stratified in com- 
position is unstable if dfj,/dz > 0, as shown in panel (b) in 
Figure 0] Because the diffusion of ions along magnetic field 
lines is inefficient, the mean molecular weight of a fluid ele- 
ment that is displaced with £ z > (£ 2 < 0) is lower (higher) 
than the surrounding medium. This results in a displaced fluid 
element with a density which is lower (higher) than the den- 
sity of the surrounding medium, which will thus rise (sink). 

If the temperature and the composition gradients are both 
positive (negative) then these arguments act in consonance 
and lead to the conclusion that the plasma is unstable (stable). 
On the other hand, if the temperature and the composition 
gradients have different signs it follows that if dT/dz > 
(d/i/dz > 0) is steep enough then the expansion induced in 
an upwardly displaced fluid element can offset the stabilizing 
effects of din fi/dz < (d In T/dz < 0) and the plasma will 
be unstable, giving rise to an HPBI. 

The arguments outlined here can also be derived from the 
equation of motions for the Lagrangian displacement £ z . The 
buoyancy force per unit volume on a vertically displaced fluid 
element produces an acceleration given by 



£6 
dt 2 



6 P 
-9— ■ 
P 



(99) 



and thus, according to Equation 



dt 2 



(100) 



This leads to an instability if iV|. > 0, in agreement with the 
results of Section |4] 



6.2. Diffusive Heat- and Particle-Flux Driven Buoyancy 
Instability 

The effects of ion-diffusion are not negligible for the modes 
for which the associated timescales satisfy t^ 1 > r^ 1 ~ 

S> Wdyn- To leading order, the changes in the tempera- 
ture and the mean molecular weight in a fluid element which 
is vertically displaced by £ z in a background magnetic field 
with b z = 1 are given by 



5T 



dT 
dz 



at* = ~riz 

dz 



(101) 



Note that because the dominant terms in Equation ([TBI are 
both proportional to D, the fractional change in the mean 
molecular weight is independent of the value of the diffusion 
coefficient. The fractional change in density is thus 



dp 
P 



d In T dln/j, 



dz 



dz 



D-HPBI. (102) 



The role played by the background temperature gradient is 
identical to the one discussed in the absence of ion-diffusion, 
and this situation is shown for the sake of clarity in panel (c) in 
FigureS] However, the contribution from the mean molecular 
weight gradient is now the opposite. This difference in sign is 
due to the manifest role played by ion-diffusion as an effective 
agent to tap into the free energy available in the background 
particle flux needed to maintain the composition gradient. 

In order to understand the role played by ion-diffusion let us 
focus on a background with constant temperature, dT/dz = 
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0. The term proportional to the composition gradient is re- 
sponsible for the onset of the D-HPBI when dfi/dz < 0. This 
is illustrated in panel (d) in Figure |4] which shows that the 
mean molecular weight of a fluid element that is displaced 
with £ z > (£z < 0) will decrease (increase) because of the 
decrease (increase) in the flux of particles out of (into) it due 
to the diverge (convergence) of field lines that results from 
the displacement of the fluid parcel. This causes the den- 
sity of the fluid element to decrease (increase) with respect 
to the surrounding medium, leading to a runaway process. A 
background temperature gradient with dT/dz > will rein- 
force this process leading to enhanced buoyancy. In general, 
when both the temperature and the composition gradients are 
non-zero, the condition for D-HPBI is din T/dz < dlnfx/dz. 
This is reflected in the equation of motion for the Lagrangian 
displacement 



ft* 
dt 2 



(103) 



which has exponentially growing solutions when -ZV|y M > 0, 



as we have seen before. 



6.3. Magneto-Thermo-Compositional Instability 

Let us now consider a background horizontal magnetic field 
with b x = 1 and modes for which conduction is faster than 
any other timescale. When a fluid element is displaced by £ 2 
from its equilibrium position, the fluctuations in temperature 
and composition are given by 



AT dT t 

dz 



dz 



(104) 



It is important to note that this relationship between the frac- 
tional change in the mean molecular weight and the La- 
grangian displacement holds regardless of whether the dy- 
namical timescale is fast or slow with respect to the viscous 
and diffusion timescales. The only difference is that in the 
former case the leading order terms in Equation (ITol that lead 
to 5/j. = — (d/i/dz)£ z are independent of the coefficient D, 
while in the opposite limit, the dominant terms are those pro- 
portional to D. In either case, the relative change in density 
becomes 



5p 



d In T dhifj, 



dz 



dz 



6 



MTCI. 



(105) 



In this case, the term proportional to the temperature gradi- 
ent is responsible for the onset of the MTI when dT/dz < 
in a homogeneous medium, i.e., dfi/dz = 0. In panel (e) in 
Figure|4]we consider the situation where r^ 1 S> Wdyn- Under 
this condition, a fluid element that is displaced with £ z > 
(£ 2 < 0) is effectively heated (cooled) by conduction along 
the magnetic field lines which have been distorted by the dis- 
placement of the fluid element. This causes the fluid element 
to expand (contract), and thus attain a density which is lower 
(higher) than the surrounding medium. This leads to the run- 
away process known as the MTI. The term proportional to the 
composition gradient is responsible for the onset of the MTCI 
when dfi/dz > in an isothermal environment. This is illus- 
trated in panel (f) in Figure|4] which shows that because of the 
effective diffusion of ions along distorted magnetic field lines, 
fluid elements that have been displaced upward (downward) 
£ z > (£, < 0) maintain the mean molecular weight cor- 
responding to the equilibrium value at their original position. 
This implies that the fluid element is immersed in a medium 
that is relatively denser (lighter) and it will thus rise (sink). 



When both gradients are non-zero their relative magnitudes 
determine wheth er th e plasma is buoyantly unstable accord- 
ing to Equation ( 1105b . This can also be seen in the equation 
of motion 



dt 2 



(106) 



which has exponentially growing solutions when N%, , < 0. 



6.4. Overstable Modes 



The physics driving overstable modes, which are present 
when b z = 1, is more subtle, but some insight can be gained 
by analyzing directly the equations of motion for the corre- 
sponding Lagrangian displacements. Using Equations d57| i 
and d78t together with the equations for the perturbations 
(IT3]>-(fT6]l we find, to first order, 



ft* 
dt 2 

dt 2 



N 2 



N 



dt 1 



(r c 1 > uj dyn ) : 



(107) 



N 2 (k 2 ± /k 2 ) dt 



,( Wdyn »r c - 1 ). (108) 



The physics behind the first terms on the right-hand side of 
each of these equations is readily recognized. The term pro- 
portional to Nj, is responsible for the HPBI if Nj, > 0, 
while the term proportional to N 2 is responsible for Brunt- 
Viasala oscillations, i.e., stable g-modes, if N 2 > 0. In a 
medium where 7V|, < and N 2 > but N 2 + 7V|, < 0, 
the anisotropic heat flow along magnetic field lines tends to 
monotonically increase the restoring force acting on the os- 
cillating fluid parcel and leads to overstability. 

7. ASTROPHYSICAL IMPLICATIONS 

Throughout this paper, we have studied the various insta- 
bilities that can be present in the parameter space spanned by 
(dln/x/dlnP, dlnT/dlnP) without imposing restrictions 
on the relative values of the gradients involved. We can now 
frame our results in the context provided by observations and 
theoretical models addressing the temperature and composi- 
tion structure of galaxy clusters. 

7.1. Temperature and Composition Profiles in the ICM from 
Models and Observations 

Modern X-ray observatories, such as Chandra and XMM- 
Newton, have made it possible to obtain the temperature 
profiles of a large number of galaxy clusters, see, e.g., 
IVikhlinin et all (120061) : ITeccardi & Molendil (120081) . This has 
enabled to show that, quite generically, the gas temperature 
increases with radius, reaching values of 5 - 10 keV at 100 
kpc, and decreases toward the outer cluster regions. In several 
cases the temperature changes by a factor of two or three on 
distances comparable to the cluster radius. 

In spite of the small mass ratio between hydrogen and 
helium, because helium is relatively abundant in the pri- 
mordial gas from which clusters form, its sedimentation 
over the life-time of the cluster has the potential to produce 
important gradients in the mean molecular weight profile. 
Obtaining observational evidence to quantify the helium 
abundances that would result form this sedimentation is very 
hard because helium is completely ionized at the charac- 
teristic temperatures of typical galaxy clusters. However, 
this information is key in order to derive physical cluster 



12 



Pessah & Chakraborty 



inner ICM intermediate ICM 



outer ICM 



0.9 



0.8 



0.7 - 



n(r\ 

- Pi' I X 

* X. 

■ / x 
* ■ 

* 


* 

\ — 

\ 

\ 

\ 

t« * 
V \ 
x * 


* 

* 

* 

* 

* 

* 

* 

* 


■ X * 
\ * 
V * — 
\ \ 
X \ 
X. % 

V 


1 1 1 1 1 1 1 1 1 


- " 



10 



- 8 



- 6 



0.01 



0.1 



1.0 



■[Mpc] 



Figure 5. Schematic representation of the mean molecular weight (red solid 
line) and temperature (black dashed li ne) profiles of a repre sentative galaxy 
cluster as suggested by observations (Vikhlinin et al. 2006) and theoretical 
models (Bulbul et al. 201 1). The regions denoted by "inner", "intermediate", 
and "outer" ICM (delimited by dotted gray lines) correspond to three different 
quarters in the (d In fi/d In P, d In T/dlnP) plane (Figures [2] [5] and panel 
(a) in Figure[6}. The mean molecular weight for a homogeneous cluster with 
primordial abundance is fj, ~ 0.6. 



properties, s uch as gas mass, t otal mass and gas mass fraction 
(see, e.g.. IQin & Wul I2000i) , and even cluster distances 
dMarkevitchll2007l) from X-ray observations. This has high- 
lighted the need to understand the efficien cy of this process 
on th eo retical and observati o nal gro u nds ([F abian & Pringlel 
19771: lAbramopoulos et all 119811: iG ilfanov & Sv unvaeyl 
IQin & Wul l200q IChuzhov & Nusseri 12001 



197" 
1984 

Chuzhov & Loeb 2004; Tamura et al.l 120041: lEttori & Fabianl 
20061: IPeng & N agai 2009J). Some of these estimates predict 
an overabundance by up to a factor of a few with respect to 
the primordial value /i ~ 0.6. 

Since the realization that the stability properties of a weakly 
magnetized medium depends critically on its thermal struc- 
ture, a large number of works have been devoted to study 
the HBI and the MTI. In spite of the fact that both instabil- 
ities grow on a dynamical timescale, the onset of the most 
relevant modes depends explicitly on the local values of the 
thermal gradients. Throughout this study, we have seen that 
if the composition of the plasma is not homogeneous, it is 
the combination of both the thermal and the mean molecu- 
lar weight gradients that decides whether the plasma is sta- 
ble or not. In order to understand how important the con- 
tributions from each of these gradients can be, we consider 
the schematic representation of the temperature and the mean 
molecular weightprofiles of the representative galaxy cluster 
shown in Figure [5] The temperature p rofile shown resemble s 
the results obtained by observations (Vik hlinin et al.l 120061) . 
whereas the mean molecular weight pro file is akin to the he - 
lium sedimentation models discusses in IBulbul et alj (1201 II) . 
which are based on analytical models for the phys ical proper- 
ties of the ICM introduced in IBulbul etaLl (120 101) . Although 
very crude, these representative profiles allow us to provide 
an estimate for the values of the logarithmic gradients for 
the temperature, ( VT) /T, and the mean molecular weight, 
(V p)/ p, which play a role in determining the stability of the 
dilute ICM. If the peak in the temperature profile occurs at a 
larger radius than the peak in the mean molecular weight pro- 
file, as shown in Figure [5] then there are three distinct regions 
defined by the signs of the temperature and the composition 
gradients. Each of these regions of the ICM will have asso- 



ciated different characteristic values, which can be estimated 
according to 

L T ii L p 



VT 



Here, AT = T out — Tj n stands for the difference between 
two values across the characteristic scale L = r out — rj n , and 
T = (T out +T n )/2 is the associated mean value, with similar 
definitions for the mean molecular weight. Using the infor- 
mation available in Figure [5] we estimate thes e cha racteristic 
values for the different ICM regions in Section 1731 

7.2. Applicability of Approximations in the ICM 

Before addressing the stability of the different regions of 
the ICM, we comment on two of the approximations that we 
have made on the geometry and strength of the magnetic field, 
which have allowed us to gain physical insight while keeping 
the problem tractable. 

We have studied in detail two special cases for the orien- 
tation of the background magnetic field, viz., either parallel, 
B = Bz, or perpendicular, B = Bx, to the gravitational 
field g. These two configurations, which have received a lot 
of attention in the related literature, have the advantage of not 
only simplifying the mathematics involved but also exposing 
in a clean way the physics driving the HBI, the MTI, as well 
as the generalizations that result from including composition 
gradients. While beyond the scope of this paper, accounting 
for more general geometries is clearly necessary in order to 
describe more realistic situations. 

An important simplification in our study is the assump- 
tion that the magnetic field is so weak that the Alfven fre- 
quency is much smaller than any other inverse timescale in- 
volved. It is worth mentioning that this is the regime ex- 
plored by a number of numerical studies add ressing both fun- 
damental aspects of the MTI and the HBI ( Parrish & Stonel 
12001 120071: IParrish & Ouataerd 120081 Latter & Kunzll2012D . 
as well as the implications that the se instabilities have for 
the long-term evolution o f the ICM dBogdanovic et all I2009t 



Ruszkowski & Ohl l2010t IParrish et al.l 120081 120091 



2010 



20121: IMcCourt et al.ll201U I2012t IKunz et alJI2012l) . For the 
modes for which magnetic tension cannot be neglected, the 
Alfven timescale can become comparable or even faster than 
the dynamical and the conduction timescales and neglecting 
wa in the dispersion relation (l33l l is not a good approximation. 

For the modes for which magnetic tension is important, the 
explicit dependence on wa might introduce new stability cri- 
teria which have not been captured by our analysis. Further- 
more, magn etic tension could affect the growth rates of the 
instabilities. Kunz (2011) provides a summary of the stabiliz- 
ing effects provided by magnetic tension on HBI- and MTI- 
unstable modes in a homogeneous medium. The main phys- 
ical effect introduced by a non-zero Alfven frequency is to 
provide a cut-off for the growth of unstable modes at parallel 
wavenumbers such that fc|| va is comparable to the growth rate 
of the most unstable modes. For the magnetic field geometries 
that we analyzed, this must also be the case even in the pres- 
ence of a composition stratification. The reason for this is that, 
when either b x = 1 or b z = 1, all the contributions introduced 
by the mean molecular weight gradient appear in the form 
d In T/dz ± d In p jdz. Thus, while the growth rates and range 
of unstable modes are affected because of the changes in the 
background composition, the effects of non-negligible mag- 
netic tension on these modes, i.e., the existence of a cut-off 
parallel wavenumber, must be similar to what has been found 
in the case of a homogeneous medium. Since the plasma j3 



THE STABILITY OF WEAKLY COLLISIONAL PLASMAS 1 3 



Table 1 

Representative Parameter Values for Different ICM Regions 



ICM Region 


VT 




p 




Local Modes in Fluid Model 

Kl /2 < fey < K~ 1/2 


Fast Conduction and < uj dyn 
1/3 < fc|| < (Z?^) 1 / 2 


Slow Conduction and oja < t c 1 

OAifiKn)- 1 / 2 < fey < 1/3 


Outer-ICM 


< 


< 


10 4 


10 1 


0.3 < fey < 3 


1/3 < fcj < 30 


0.003 < fey < 1/3 


Interm.-ICM 


> 


< 


10 3 


10 2 


0.1 < fej| < 10 


1/3 < fey < 3 


0.03 < fey < 1/3 


Inner-ICM 


> 


> 


10 2 


10 3 


0.03 < fey < 30 


1/3 < fey < 0.3 


0.3 < fey < 1/3 



Note. — The various ICM regions are defined in Figure [5] For convenience, in this table we have defined the dimensionless parallel wavenumber, 
fe|| = fc||(A mfp /T) 1/2 . 



in galaxy clusters varies with radius, whether neglecting mag- 
netic tension is a sensitive approximation for a local stability 
analysis or not, depends on the conditions present in the re- 
gion of the ICM under consideration. We address this issue in 
detail below. 

7.3. Stability of ICM regions 

The analysis of Figures|2]and[3]allows us to understand the 
implications that a mean molecular weight gradient can have 
for the various regions of a representative galaxy cluster as 
depicted in Figure These regions correspond to different 
quadrants in the (d In old In P, d In T/d In P) plane as shown 
in panel (a) of Figure [6] which we have denoted as inner, in- 
termediate, and outer ICM. In what follows we will assume, 
as suggested by observations of galaxy clusters, that the ICM 
is buoyantly stable according to the classical stability criterion 
N 2 >0 (IPiffaretti et al.ll20^ ICavagnolo et al.ll20"o1l . 

7.3.1. Outer ICM 

In this region the temperature and the mean molecular 
weight gradients are both negative. Because the criteria for 
stability are different whether conduction is fast or slow com- 
pared to the dynamical timescale, we consider these two cases 
separately. 

Fast conduction. If the magnetic field lines are parallel to 
the gravitational field, i.e., b z = 1, gravity modes can be- 
come overstable if d In T/d In P > (7 — l)/(27), panel (a) 
in Figure|2] A gradient in the mean molecular weight alone is 
unable to stabilize these modes and can drive unstable modes 
driven by diffusion if |V/x|//i > \VT\/T, panels (b) and (c) 
in Figure [2] For magnetic field configurations that are per- 
pendicular to the gravitational field, i.e., b x = 1, this region is 
unstable to the MTCI provided that I V/^ I //i < |VT|/T,panel 
(d) in Figure 12 This means that the outskirts of galaxy clus- 
ters that would be considered prone to the MTI (if they were 
homogeneous) would remain stable if the gradient in mean 
molecular weight is steep enough. This is not the case for the 
particular profiles shown in Figure [5]but this does not imply 
that this is not the case in general. 

Slow conduction. In the absence of ion-diffusion, gravity 
modes can become overstable if bz = 1 . These modes cannot 
be stabilized by means of a gradient in the mean molecular 
weight alone, panel (a) in Figure [3] Furthermore, when ion- 
diffusion is efficient, it can drive unstable modes if | V/i|//i > 
|VT|/T, panel (c) in Figure|3] For b x = 1 there can be unsta- 
ble modes driven by conduction if |V/i|//i < | VT|/T, panel 
(d) in Figure [3] 

We can provide a crude estimate of the impact that a com- 
position gradient would have on the growth rates of the vari- 
ous instabilities discussed by estimating the temperature and 
the composition gradients shown in Figure [5] For the inner 



ICM region, the characteristic scale is L ~ 0.8 Mpc, while 
AT ~ 4 keV, T ~ 7 keV, A/i ~ 0.1, and jl ~ 0.65. Thus 
according to Equations (|109l l, the characteristic logarithmic 
gradients in this inner region are given by 



VT 



-O^Mpc" 1 , — 



/' 



-0.2 Mpc" 1 .(110) 



These order-of-magnitude estimates, based on the represen- 
tative values drawn from Figure show that the instabilities 
with growth rates a 2 cx — \n(T/uX such as the generalization 
of the MTI, Equations (|59]) and (1681 . will be 15% slower with 
respect to the homogeneous case. 

Regarding the validity of our assumptions of weak mag- 
netic fields, panel (b) of Figure [6] shows that magnetic ten- 
sion is unimportant for the range of modes of interest, i.e., 

K}J 2 < ku (A^pii") 1 / 2 < Kn . Thus, our approximation 
of setting oja ~ is fully justified in this region. 

7.3.2. Intermediate ICM 

In this region the temperature and the mean molecular 
weight gradients have different signs, Q with V/i < and 
VT > according to the profiles shown in Figure 

When b z = 1, this region is unstable due to the HPBI, 
which grows on the dynamical timescale if |V/i|//i < 
I VT|/T. This instability can be prevented if the mean molec- 
ular weight is steep enough, panel (a) in Figure [2] When 
ion-diffusion is efficient this region is unstable due to the 
D-HPBI, panel (c) in Figure [2] This magnetic field config- 
uration is also prone to unstable modes for which conduc- 
tion is slow. In this case there are unstable modes driven by 
conduction if D = and |V//|//x < |VT|/T, while there 
are unstable modes driven by ion-diffusion if D =^ and 
|V/x|//i > |VT|/T, panels (b) and (c) in Figure |2] respec- 
tively. If b x = 1, this region is stable whether conduction is 
fast or slow compared to the dynamical timescale, panels (d) 
in Figure|2]and (e) in Figure|3] 

For this intermediate ICM region, the inspection of Figure|5] 
provides L ~ 0.15 Mpc, AT ~ 2 keV, T ~ 8 keV, Afj, ~ 0.1, 
and /2 ~ 0.75. Thus the characteristic logarithmic gradients 
are given by 



VT 



1.6 Mpc" 1 , — 



intcrm. 



11 



1 Mpc" 1 .(111) 



intcrm. 



Therefore, the instabilities with growth rates for which a 2 cx 
ln(T^), such as the generalization of the HBI in the absence 

7 Note that if the peak in fx were to occur beyond the radius at which 
the temperature is maximum, then this region would be characterized by 
dT/dz < and dujdz > 0, which corresponds to the opposite quadrant in 
panel (a) of Figure|6] 
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Figure 6. Panel (a) shows a schematic representation of the correspondence between the regions of a representative galaxy cluster with the radial temperature 
and the mean molecular weight profiles as shown in Figures |3]and the plane spanned by (d In T/d\n P, d\n fi/d In P). The various vertical lines in panels 
(b)-(d) mark some relevant values of the dimensionless parallel wavenumber, fen (AmfpH) 1 / 2 . The region between the dashed lines con'esponds to the range 

of local modes for which the fluid approximation is valid, i.e., kU' 1 < fey (A m f p i?) 1//2 < K n The solid line, at fell (A m f p //) 1//2 = 1/3, distinguishes 
between modes for which conduction timescale is smaller (to the right) or larger (to the left) than the dynamical timescale. The dotted lines represent the value of 
the feii (A m f p i? ) 1//2 beyond which magnetic tension cannot be neglected. This con'esponds to fey (^mipH) 1 ^ 2 > (pKn) 1 / 2 (green) for fast conduction modes, 
and fen (A m f p -ff) 1//2 < (/3A'„)~ 1 / 2 /10 (blue) for slow conduction modes. All the numerical values associated with the various length scales and timescales 
have been drawn from Table[T] which provides estimates for the f) plasma and the Knudsen number representative of the different regions of the ICM. 



of ion-diffusion, Equation ( 154b . will grow of the order of 40% 
slower with respect to the homogeneous case. On the other 
hand, the instabilities with growth rates a oc ln(T//j,), such 
as the generalization of the HBI with ion-diffusion, Equa- 
tion ((63, will be 30% faster. 

As shown in panel (c) in Figure [6] magnetic tension is im- 
portant for some of the modes for which conduction is fast but 
not for the modes for which conduction is slow. 

7.3.3. Inner ICM 

In this region the temperature and the mean molecular 
weight gradients are both positive. We consider again the lim- 
its in which conduction is fast or slow separately. 

Fast conduction. For b z = 1, this region is unstable to the 
HPBI regardless of whether the mean molecular weight gra- 
dient is smaller or larger than the temperature gradient, panel 
(a) in Figure[2] Furthermore, if ion-diffusion is efficient it can 
also drive unstable modes, panels (b) and (c) in Figure [2] In 
a homogeneous medium with b x = 1, this inner region is sta- 



ble against the MTI. However, there can be unstable MTCI- 
modes if |V>u|//i. > |VT|/T, panel (d) in Figure|2] In the ho- 
mogeneous case when VT > 0, the HBI tends to re-orient the 
magnetic field in the radial direction, which results in a field 
configuration which is stable against the MTI, i.e., b x ~ 1. 
When ion-diffusion is not efficient, this core insulation could 
be alleviated by the MTCI if the mean molecular weight gra- 
dient is steep enough, i.e., (V/x)//i > (VT)/T. 

Slow conduction. This region can be subject to instabili- 
ties driven by both heat conduction and ion-diffusion. When 
b z — 1 and ion-diffusion is inefficient, there are unstable 
modes driven by heat conduction regardless of the relative 
magnitude of the temperature and the mean molecular weight 
profiles, panel (b) in Figure [3] whereas there are unstable 
modes driven by diffusion if |V/i|//i > |VT|/T, panel (c) 
in Figure [3] In the case with b x = 1, there can be unsta- 
ble modes driven by heat conduction if |V/i|//i > \VT\/T, 
panel (e) in Figure [5] 

According to Figure [5] the inner ICM region is charac- 
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terized by L ~ 0.05 Mpc, AT 
A/i ~ 0.1, and p, ~ 0.8, and thus 



10 Mpc" 



V/i 

I" 



3 keV, T ~ 5.5 keV, 



2. 5 Mpc" 1 . (112) 



Therefore, the instabilities with growth rates for which cr 2 oc 
ln(T/i) will grow of the order of 10% faster with respect to the 
homogeneous case, while the instabilities with growth rates 
cr 2 oc ln(T//i) will be 15% slower. 

The analysis of Figure|6]shows that the range of modes for 
which it is sensible to carry out a local mode analysis within 
the fluid model embodied in Equations (d}-© increases as 
the inverse Knudsen number increases toward smaller radii. 
However, due to the increase in the strength of the background 
magnetic field, the range of modes for which it is sensible to 
neglect magnetic tension, decreases. Because the relatively 
low values of j3 in the inner core region, magnetic tension is 
important for all the modes of interest. Therefore, our conclu- 
sions for this region should be considered with caution. 

7.3.4. Summary and Outlook 

During the last few years, there has been substantial 
numerical work for understanding the long-term evolu- 
tion of the MTI and the HBI and their implications for 
the gas dynamics in the ICM permeating g a laxy c lusters 
dBogdanovic et alj l2009t IParrish et all 120081 120091 l20T0i 
120121; IMcCourtetal.1120111 120121: IKunz et alj 120121) . All 
of this work has been done under the assumption that 
the ICM is homogeneous and thus the temperature gradi- 
ent provides the only source of energy to feed instabili- 
ties. Even though it is hard to quantify concentration gra- 
dients from observations, some heavy element sedimentation 
is expected dNaravan & M edvedev 200 11 IChuzhov & Nusserl 
12001 IChuzhov & Loebll2004t lEttori & Fabianl2006l) . Indeed, 
current theoretical models suggest that helium sedimenta- 
tion can significantly alter the composition profile throughout 
the cluster and give rise to mean molecular weight gradients 
which are comparable in magnitude to the temperature gradi- 
ents, with |VT|/T ~ |V//|//x (see, e.g.. IBulbul eTaLll20Tll 
and Figure [5J. 

This work discusses for the first time the effects that com- 
position gradients can have for the stability of a weakly colli- 
sional magnetized medium which is stratified in both temper- 
ature and composition. We have found that, depending on the 
wavelength of the modes under consideration, the plasma can 
be subject to a wide variety of unstable modes. These include: 



the generalization of the MTI ( IBalbusll200Tb : 

dln(T/ii) k 2 x 



a K,-g- 



k y 



dz k 2 ' 

and the generalization of the HBI (jpuataert 2 0081) in 

• the slow ion-diffusion limit: 

2 dln(7» k 2 

• the fast ion-diffusion limit: 

g dhi(T j \x) 
tJ 1 dz 



(113) 



(115) 



We have also found th e generalization of the over stable grav- 
ity modes discussed in lBalbus & Reynolds! ( 120101) . see Equa- 
tions d57| i and ( 17 81 . as well as other new modes which are 
driven by conduction and diffusion. 

This study constitutes a first step toward the long-sought 
goal of understanding in a self-consistent way the effects of 
magnetic turbulence on the diffusion of heavy elements and 
its consequences for the observational signatures and long- 
term evolution of galaxy clusters. This will only be possible 
through numerical studies involving realistic models for the 
microphysics of weakly collisional, multi-component plas- 
mas. Addressing this problem will require to sort out sev- 
eral details, including how to properly handle plasma micro - 
instabilities, many of which are still the subject of active re- 
search in homogeneous settings. 
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APPENDIX 

A. ION DIFFUSION IN A BINARY MIXTURE 

The ratio between the timescales associated with viscous 
and diffusion processes is t~ 

X I t a 1 = 3v\\/D. Here, the co- 
efficient i/ii = v 2 h /(2is^) denotes the kinematic Braginskii 

viscosity associated with a binary mixture of ions, where i/?f 
is an effective ion-ion collision rate which can be estimated 
as follows. The Braginskii viscosity for a single species of 
ions is ?7o = pv 2 h /(2z^), where va is the ion-ion collision 
frequency 



4y^ 



m\ /2 {k B Tf/ 2 



In A, 



(Al) 



Here n^, and qi are the number density, the mass, and the 
charge of the ion respectively; and In A refers to the corre- 
sponding Coulomb logarithm. In a binary mixture, the ef- 
fective viscosity coefficient (ignoring the contribution from 
electrons) is given by 



n ix k B T 



n l2 k B T 
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where the collision frequency between ions of species i and j 
is 



m^ksTf/ 2 



In A 
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and rriij = mimj/(m 



rrij) is the reduced ion mass. We 



(1 14) can thus define the effective ion-ion collision frequency as 
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The coefficient governing the diffusion of species 2 
into species 1 (e.g., helium into hydrogen) is given by 



16 
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Figure 7. Ratio of the inverse timescales associated with viscous and diffusion processes for a binary mixture of H and He as a function of the composition c 
(left) and mean molecular weight fi (right). The vertical dashed line indicates the values for which the mixture has primordial composition, i.e., c ~ 0.25, which 
corresponds to fi ~ 0.6. 
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Therefore the ratio t v 1 /r d 1 = ivnjD is given by 



(2-c)(8-5c) 
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which is shown in Figure|7]as a function of the concentration c 
(left panel) and as a function of the mean molecular weight /i 
(right panel). In the inner regions of the ICM, where c ~ 0.6 
or /i ~ 0.8 (see Figure 0, viscous and diffusion processes 
take place on comparable timescales, i.e., r^ 1 /t^ 1 ~ 3, 

while t~ 1 /t7 1 ~ 9 for a primordial mixture of helium and 
hydrogen, i.e., c ~ 0.25 or /i ~ 0.6. It should be kept in mind 
that in carrying out this calculation we have assumed that all 
the ratios between the Coulomb logarithms are of order unity 
and the expressions for the transport coefficients along the 
magnetic field lines are identical to the ones that are valid in 
the absence of the magnetic field. 
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